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Determining the cellular composition of specific brain regions is crucial to our 
understanding of the function of neurobiological systems. It is therefore useful to identify 
the extent to which different methods agree when estimating the same properties of 
brain circuitry. In this study, we estimated the number of neuronal and non-neuronal cells 
in the primary visual cortex (area 17 or V1 ) of both hemispheres from a single chimpanzee. 
Specifically, we processed samples distributed across V1 of the right hemisphere after 
cortex was flattened into a sheet using two variations of the isotropic fractionator cell 
and neuron counting method. We processed the left hemisphere as serial brain slices 
for stereological investigation. The goal of this study was to evaluate the agreement 
between these methods in the most direct manner possible by comparing estimates 
of cell density across one brain region of interest in a single individual. In our hands, 
these methods produced similar estimates of the total cellular population (approximately 
1 billion) as well as the number of neurons (approximately 675 million) in chimpanzee 
V1, providing evidence that both techniques estimate the same parameters of interest. In 
addition, our results indicate the strengths of each distinct tissue preparation procedure, 
highlighting the importance of attention to anatomical detail. In summary, we found that 
the isotropic fractionator and the stereological optical fractionator produced concordant 
estimates of the cellular composition of V1, and that this result supports the conclusion 
that chimpanzees conform to the primate pattern of exceptionally high packing density 
in V1. Ultimately, our data suggest that investigators can optimize their experimental 
approach by using any of these counting methods to obtain reliable cell and neuron counts. 
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INTRODUCTION 

Determining the cellular composition of neurobiological tissue is 
critical to our understanding of the structure and function of the 
brain. Differences in the organization of neurobiological tissue, 
whether it is the arrangement of cells or the relative amount of 
an organism's nervous system dedicated to a particular sensory 
modality, are important components of the neurobiological cir- 
cuitry underlying behavior (Sur et al, 1980; Catania et al., 1993; 
Krubitzer et al, 1995). 

Scientists have employed a variety of techniques with varying 
degrees of success to more fully characterize the microstructure 
and circuitry of specific brain areas and systems. An early exam- 
ple of such a procedure is "direct enumeration," which consists 
of homogenizing brain tissue and determining the density of cell 
nuclei in re-suspended fluid samples (Nurnberger and Gordon, 
1957; Zamenhof, 1976). A major critique of this approach was 
concern over the potential loss of cells or nuclei during the 
mechanical dissociation of tissue (Clarke and Oppenheim, 1995). 
Alternative procedures were developed to analyze specific param- 
eters of sectioned tissue samples, and since the 1960's these 
stereological techniques have become the gold-standard for quan- 
titative investigations of biological tissues (Glaser and Glaser, 
2000). However, stereological methods have also been critiqued 
for inattention to tissue processing deformations that may lead 



to inaccurate estimates (von Bartheld, 2002). Traditionally, for 
example, guard zones were used to avoid tissue containing 'lost 
caps" that resulted from cells being destroyed by direct contact 
with the knife blade (Baryshnikova et al., 2006). However, recent 
work suggests the absence of "lost caps" in fixed tissues cut on a 
freezing microtome, but reports that the distribution of particles 
along the z-axis (from one cut edge to the other) is non-uniform, 
concluding that the traditional use of guard zones may intro- 
duce bias into the estimate for these tissues (Carlo and Stevens, 
2011). Indeed, tissue deformations that introduce bias into stere- 
ological investigations have been reported for a number of tissue 
processing protocols and constitute a major cause for concern 
(von Bartheld, 2012). Thus, the extent to which these distinct 
methodologies produce accurate estimates of specific biological 
parameters of interest remains controversial (Bahney and von 
Bartheld, 2014). Accordingly, whether these methods produce 
comparable results has also been of interest for several decades as 
agreement between methods provides support for the conclusion 
that we are, in fact, producing reasonable estimates of a specific 
biological parameter of interest, and forms the motivation for our 
work. An additional challenge to this type of work is the abil- 
ity to accurately identify the region of interest for quantitative 
analysis. Immunohistochemical preparations have enhanced our 
ability to identify discrete brain areas by providing information 
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on an increased number of the molecular properties of the brain 
(Sherwood et al., 2007a,b; Raghanti et al, 2008; Takahata et al., 
2009; Balaram et al., 2013). Indeed, stereological analysis has 
seen widespread success in part because of the ability to employ 
immunohistochemical protocols to tag neurons (Mullen et al., 
1992) as well as histological procedures to visualize cells, includ- 
ing the Nissl stain (Miller et al, 2012), and neurotransmitter 
biomarkers, such as acetylcholinesterase (Miller et al., 2013), in 
adjacent brain slices. Not only does this increased number of 
staining procedures give researchers more to investigate with the 
techniques at their disposal, it also provides an enriched view 
of the neurobiology under investigation and facilitates the rapid 
and reliable identification of brain areas. Importantly, recent years 
have witnessed the application of some of these newer technolo- 
gies to older procedures. For example, the isotropic fractionator 
(Herculano-Houzel and Lent, 2005) is the product of combining 
"direct enumeration" with immunohistochemistry to investigate, 
for example, the number of neuronal and non-neuronal cells in 
a tissue sample. In addition, results from manual counts using 
the isotropic fractionator method (Herculano-Houzel and Lent, 
2005) have been replicated using modern flow cytometry to 
more rapidly quantify the cellular composition of neurobiolog- 
ical samples, termed the flow fractionator (Collins et al., 2010; 
Young et al, 2012). A principal reason these procedures have seen 
widespread use in recent years is that each has specific advan- 
tages. For instance, while stereological analysis of brain slices 
provides a wealth of anatomical detail, the preparation of his- 
tological sections is time-consuming, and adequately sampling 
brain tissue can be challenging. Conversely, the isotropic frac- 
tionator can more rapidly produce reliable estimates of some of 
the same parameters, such as the cellular composition of a brain 
sample, via homogenization of samples (Herculano-Houzel and 
Lent, 2005). Recently, the use of flow cytometry to automati- 
cally count samples prepared using the isotropic fractionator was 
demonstrated to further decrease the amount of time necessary 
to estimate cellular populations, and to provide results less vari- 
able than those obtained from hand counts using a microscope, 
although with an increased equipment cost (Collins et al., 2010; 
Young et al., 2012). Furthermore, although stereological studies 
can employ a range of histological markers to identify individual 
regions of interest, the isotropic fractionator and flow fractiona- 
tor can each be used for samples taken from cortex that has been 
separated from the white matter and flattened into a sheet, in 
which specific brain areas, such as the primary visual cortex (VI), 
are readily identifiable during dissection. Despite the theoretical 
ability of these procedures to estimate the same parameters, no 
study has to-date directly compared results from the use of the 
isotropic fractionator, the flow fractionator and the stereological 
optical fractionator in a single brain region. 

The present study is therefore the first to investigate the sim- 
ilarity of estimates obtained using the isotropic fractionator, the 
flow fractionator and the stereological optical fractionator from 
a specific region of the brain within a single individual. We 
determined the total number of cells by staining free-floating sus- 
pensions of nuclei for DNA with 4',6-diamidino-2-phenylindole 
(DAPI), and the percentage of neuronal cells by staining suspen- 
sions of nuclei with both DAPI and an antibody to the neuronal 



nuclei marker NeuN (Mullen et al., 1992) in samples from flat- 
tened VI of the right hemisphere using the isotropic fractionator 
(Herculano-Houzel and Lent, 2005; Collins et al, 2010; Young 
et al., 2012). We then determined the total number of cells and 
the percentage of neuronal cells in the same samples with the 
flow fractionator (Collins et al., 2010; Young et al, 2012). Finally, 
we recorded the location of cells along the z-axis to determine 
the presence and extent of possible non-uniform tissue shrinkage 
that could bias our stereological estimates, and calculated the total 
number of cells and neurons using the stereological optical frac- 
tionator (Schmitz and Hof, 2005) in adjacent series of brain slices 
from the left hemisphere stained for Nissl substance or NeuN, 
respectively. Determining the agreement in results between these 
methods is important because each of these methods are in cur- 
rent use and the reliability of our conclusions and inferences 
based upon such data will depend upon how they compare. In 
addition, if estimates from these methods agree, this will pro- 
vide investigators with a greater degree of experimental flexibility 
and increase the data to which they have access when interpret- 
ing their results. In short, determining the transparency of results 
from these techniques facilitates the integration of a larger set of 
data in the development of formal relationships between the cel- 
lular composition and behavioral function of brain regions across 
a broad range of taxa and in both clinical and basic scientific 
applications. 

METHODS 

SAMPLE PREPARATION 

Both hemispheres of the brain from a single, non- 
neuropathological adult female chimpanzee (Pan troglodytes; 
53 years) were obtained from the Texas Biomedical Research 
Institute's Southwest National Primate Research Center. The 
brain was collected after post-mortem flushing with 0.1% phos- 
phate buffered saline (PBS), and was shipped overnight in 0.1% 
PBS. Upon arrival, the cerebral cortex was separated from the 
subcortical structures, bisected, and then cleaned of vasculature 
and pia. In the right hemisphere, the sulci were opened in order 
to flatten the cortical sheet, which was then fixed in 2% PFA in 
PBS for 72 h at 4°C. The left hemisphere was cut into blocks and 
fixed in 4% PFA in PBS for 24 h at 4°C, followed by immersion 
in 30% sucrose overnight prior to sectioning at 50 \im on a 
sledge microtome (American Optical, model #860, clearance 
angle =19°) while the tissue was held below 0°C. 

Details of the tissue processing methods employed in flattening 
cortex and preparing samples for the isotropic fractionator pro- 
tocol appear in prior publications (Herculano-Houzel and Lent, 
2005; Collins et al, 2010; Young et al, 2012). Briefly, the flattened 
(right) cortical hemisphere was photographed and placed onto a 
light box. Primary visual cortex (VI or Brodmann's area 17) in 
flattened cortex is easily identifiable when backlit due to its dense 
myelination and caudal position. Identifiable cortical areas (VI) 
were drawn on the printed image of the flattened cortex and VI 
was separated from the rest of cortex and subdivided into pieces 
using a scalpel. Each piece was photographed, weighed, and doc- 
umented as a single sample. Photographs were taken in order to 
measure the surface area of each piece using ImageJ, which were 
combined with the average cortical thickness of VI based upon 
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measurements (mean depth = 1.8 mm, n = 42) in serial sections 
from the left hemisphere to determine the volume of VI. Each 
sample from the flattened cortical tissue (n = 61) was homoge- 
nized using a glass Tenbroeck tissue grinder (Fisher Scientific) in 
a dissociation solution of sodium citrate, triton-X 100 and dis- 
tilled water. The total suspension volumes were determined based 
upon sample weight, and ranged between 3 and 5 ml. The result- 
ing solution contains free-floating nuclei, which were stained with 
DAPI, which indiscriminately stains cellular nuclei, and an anti- 
body to Neuronal Nuclei (NeuN, MAB377, Millipore), which 
stains most neuronal nuclei (Mullen et al., 1992). The isotropic 
fractionator and flow fractionator procedures described below 
permit the estimation of total cell number as well as the per- 
centage of neurons in the sample, assuming that each nucleus 
corresponds to one cell. 

ISOTROPIC FRACTIONATOR 

Nuclei in samples from the flattened (right) VI and processed 
using the isotropic fractionator were stained with DAPI and 
an antibody to NeuN and quantified by hand in a Neubauer 
counting chamber under a fluorescent microscope (Herculano- 
Houzel and Lent, 2005). Specifically, free-floating, DAPI-stained, 
or DAPI+/NeuN+ double-stained nuclei from the main sus- 
pension samples were counted to estimate total cell number or 
percentage of neurons, respectively (Figure 1, Table 1). Initial 
manual counts of DAPI-stained or DAPI+/NeuN+ double- 
stained nuclei were performed by a technician blind to the 
flow cytometry data, and control samples revealed excellent 
agreement between the initial and final manual counts per- 
formed by DJM on DAPI-stained nuclei (nonparametric T-test, 
P = 0.967, n = 17), and between manual and flow counts of 
DAPI+/NeuN+ double-stained nuclei (nonparametric T-test, 
P = 0.123, n = 13). All suspensions were vigorously mixed prior 
to chamber loading and counting, and each of the 61 samples 
from right VI was processed using the isotropic fractionator 
method. 

FLOW FRACTIONATOR 

Nuclei in samples from the flattened (right) VI and processed 
using the flow fractionator were stained with DAPI and an anti- 
body to NeuN and quantified by a flow cytometer (Collins et al., 
2010; Young et al., 2012). The flow fractionator also requires the 
addition of a known number of Countbright Beads (Invitrogen) 
(Young et al, 2012) (Figure 1). Specifically, free-floating, DAPI- 
stained or DAPI+/NeuN+ double-stained nuclei from the main 
suspension samples were counted in tandem with a fixed vol- 
ume of fluorescent Countbright beads (Invitrogen Inc.) using 
a Becton Dickson 5-laser LSRII flow cytometer (for details, see 
Collins et al, 2010; Young et al., 2012) (Figure 1). All samples 
were vigorously mixed before counts were determined, and all 
counts were made in duplicate to assess variability. The flow 
cytometry technician was blind to the sample attributes pre- 
pared by DJM and to the manual count data. All flow cytometry 
experiments were conducted in the Vanderbilt University Medical 
Center Flow Cytometry Shared Resource. Each of the 61 sam- 
ples from right VI was processed in duplicate using the flow 
fractionator method. 




DAPI-A 

FIGURE 1 | Summary figure of the optical fractionator, isotropic 
fractionator, and flow fractionator. Low magnification i~1x, scale bar = 
500 \im) images of coronal brain slices at the V1/V2 cortical boundary 
(indicated by the arrowhead) stained for Nissl (A) and an antibody to NeuN 
(B) showing the characteristic laminar profile of V1 that was used to guide 
stereological measurements. High magnification images (100x, scale bar = 
10 |im) of the outer stripe of Gennari in V1 in coronal brain slices stained for 
Nissl (C) and an antibody to NeuN (D), with the stereological probe 
illustrating lines of inclusion (green) and exclusion (black). High 
magnification (20x, smallest individual boxes are 50 by 50 u.m, scale bar = 
40|j,m) images of flattened cortical samples processed using the isotropic 
fractionator and quantified using the Neubauer chamber (E,F,G). (E) Shows 
DAPI-stained nuclei, (F) shows NeuN-stained nuclei and (G) is a merged 
image of (F,G). Yellow arrowheads indicate nuclei double-positive for 
fluorescent staining of DAPI and NeuN, indicating neuronal nuclei, and blue 
arrowheads indicate nuclei that stained for DAPI, but did not stain for 
NeuN, indicating non-neuronal nuclei (E,F,G). (H) Shows an example of our 
data generated using the flow fractionator and presented on a SSC-A (side 
scatter area) vs. DAPI-A (DAPI-fluorescent area) scatterplot. A polygon gate 
was used to select the nuclei based on DAPI expression. This analysis 
profile allowed us to quantify the concentration of individual nuclei (singlet, 
teal), clusters of two nuclei (doublet, red) and clusters of three nuclei 
(triplet, black) in a known volume of sample. The green dots on the left side 
of the plot indicate debris excluded by the polygon gate (black diagonal 
line). 
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Table 1 | Adult chimpanzee visual cortex total cell and neuron estimates based on three quantification techniques. 



Tissue 


Quantification 


Parameter 


Sample 


Estimate 
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Neubauer chamber 
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n = 61 


nno a on 1 /i o 
yyo,4oU, l4o 


Mat 


Flow cytometry 


All ^nllo 

All ceils 


n — o i 


i nic ccc o/io 
1 ,U 1 0,DDD,o4y 


Slice 


Optical fractionator 


All ^^llr- 

All cells 


n = 324 


yb 1 ,(Job,4bU 


Mat 


Neubauer chamber 


Neurons 


n = 61 


CC1 "7QO O 1 /I 

bb 1 , /oy,z 14 


rlat 


Flow cytometry 


Neurons 


n = 61 


cc /i tic no -i 

bb4, /zb,yo 1 


Slice 


Optical fractionator 


Neurons 


n = 316 


enc /i~7o ocn 
byb,4/o,zbU 


rlat 


Neubauer chamber 


Percent neurons 


n = 13 


62.0% 


riat 


Flow cytometry 


Percent neurons 


n — o i 


bb. I 7o 


Slice 


Optical fractionator 


Percent neurons 


n = 10 


72.4% 


Flat 


Image J 


Volume (mm 3 ) 


n = 61 


5275 


Slice 


Cavalieri (Nissl) 


Volume (mm 3 ) 


r-i 1 n 

n = iu 


4o iy 


Slice 


Cavalieri (NeuN) 


Volume (mm 3 ) 


n= 10 


4778 


Flat 


Neubauer chamber 


Cell density 


n = 61 


192,184 


Flat 


Flow cytometry 


Cell density 


n = 61 


195,551 


Slice 


Optical fractionator 


Cell density 


n = 10 


197,331 


Flat 


Neubauer chamber 


Neuron density 


n = 61 


125,641 


Flat 


Flow cytometry 


Neuron density 


n = 61 


128,005 


Slice 


Optical fractionator 


Neuron density 


n = 10 


146,008 



Summary table showing the parameter estimates obtained from flattened cortical samples (Flat) and brain slices (Slice). 



OPTICAL FRACTIONATOR 

This study also employed a traditional histological procedure to 
reveal cell bodies based upon the presence of staining for Nissl- 
substance (Nissl), as well as immunohistochemical staining for 
NeuN (Young et al., 2013), in which VI is easily identifiable 
(Figure 1), for stereological quantification. Specifically, adjacent 
l-in-10 series of sections cut from the left hemisphere on a freez- 
ing sledge microtome were processed for Nissl-substance using 
1% thionin acetate (T3387 Sigma-Aldrich) dissolved to 0.025% 
into an acetate buffer (containing sodium acetate and acetic acid 
to a pH of 5.1-5.5) after mounting on glass slides and dried 
overnight in an oven set to 37°C, or free-floating with an anti- 
body to NeuN. Sections were dehydrated in a graded series of 
alcohol (70, 95, 100%) followed by immersion in xylene prior 
to being coverslipped with Permount (SP15-500 Fisher). Neurons 
were apparent in both Nissl- and NeuN-stained materials because 
of the presence of a round nucleus. Glia were apparent in Nissl- 
stained sections because of their smaller size and darkly stained 
punctation. All stereological analyses were performed by a single 
investigator (DJM) using a computer-assisted stereology system 
consisting of a Nikon E80i light microscope with a motorized 
stage and the commercially available Stereolnvestigator software 
(MicroBrightField Bioscience, Williston, VT). The volume of each 
region of interest was estimated using the Cavalieri principle 
(frame area = 500 x 500 u,m) at 4x magnification (Table 1). An 
initial pilot study, in which sections were outlined at 2x mag- 
nification and counts were placed at 100 x magnification using 
an oil immersion objective lens under Kohler illumination when 
the top of the nucleus, in the case of neurons, or the top of 
the cell profile, in the case of glia, came into focus within the 
probe dimensions (Schmitz and Hof, 2005), was performed to 
determine stereological parameters. Specifically, the pilot study 



indicated the appropriate number of sections to be investigated 
(n = 10), probe area (40 x 40 [im), probe height (10 (Jim), dual 
vertical guard zone height (2 (Jim), frame step (1966 x 2595 (xm), 
frequency of thickness measurement (1 in 8; measured thick- 
ness: mean = 16.8 (xm, s.e.m. = 0.31 (xm) while maintaining 
an acceptably low coefficient of error [Gundersen (m = 1) CE 
< 0.025; 1st Schmitz and Hof CE < 0.024; 2nd Schmitz and Hof 
CE < 0.019] (Gundersen and Jensen, 1987; Schmitz and Hof, 
2005). 

In addition, these preliminary data indicated that the chro- 
magens used to visualize cells or nuclei achieved full penetration 
of the tissue, based upon the homogeneity of cell counts across 
the depth of probes used to investigate NeuN-stained materi- 
als, and the failure to reject the null hypothesis, in Nissl-stained 
materials, that there was no difference in the number of cells allo- 
cated to sequential counting bins along the z-axis [-F(i,8) = 0.82, 
P > 0.39, n = 10 bins for a probe depth of 10 u,m]. Additional 
analysis of the vertical (z) distribution of particles was performed 
at 100 x magnification by recording the z-axis location of cells 
(nuclei; n = 1399) within the measured depth of tissue at sam- 
pling sites (n = 207) from across VI in Nissl-stained material. 
A total of 640 sampling sites were investigated using the optical 
fractionator workflow to estimate cell and neuron number (mean 
number of sites per section = 32.4, 34.2; standard error of the 
mean (SEM) = 5.8, 4.8; NeuN, Nissl, respectively). The num- 
ber of cells or neurons was calculated by multiplying the number 
of counts by the reciprocals of the area section fraction (asf), 
serial section fraction (ssf), and the tissue section fraction (tsf) 
(Schmitz and Hof, 2005). The percent of neuronal cells in brain 
slices investigated using stereology was calculated by taking the 
number of cells estimated using NeuN-stained sections divided 
by the number of cells estimated using Nissl-stained sections. 
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STATISTICAL ANALYSIS 

The goals of this project were (1) to determine if manual and 
automated methods of counting single and double-stained nuclei 
in suspension volumes of homogenized tissue samples resulted 
in statistically indistinguishable estimates of total cell and neuron 
numbers from readily identifiable regions of interest in flattened 
cortical preparations, (2) to compare these results with estimates 
of the same parameters in the same region of interest from the 
other hemisphere in a single individual using stereological proce- 
dures on sectioned tissue, and (3) to determine the distribution 
of particles in the z-domain of Nissl-stained sections. Pursuant to 
the first goal, we determined manual counts as well as automated 
flow counts from at least two aliquots from each homogenized 
sample of flattened cortex used in this study (n = 61, Table 1). 
Next, we determined within-measures agreement by calculating 
the non-parametric T statistic on repeated measurements of each 
aliquot, testing the difference between means, such that a P-value 
smaller than 0.05 was taken to indicate a significant difference 
between the measurements (Table 2). We then used the Bland and 
Altman (1986) approach of comparing the differences between 
the measurements to the means of the measurements in our anal- 
ysis of the isotropic fractionator and the flow fractionator data 
(Figures S1-S4). In addition, we report Lin's concordance corre- 
lation, a statistic that approximates the variation of the data about 
the linear regression (Lin, 1989) (Table 2). Pursuant to the second 
goal of comparing results from samples of flattened cortex with 
results from coronal brain slices, we calculated the density of cells 
or neurons per cubic millimeter by taking the parameter estimate 
of the sample divided by the measured volume of that sample. 
For example, when using the optical fractionator, the investigator 
draws an outline of the region of interest in each of the 10 sections 
selected for analysis. The Cavalieri method then utilizes this out- 
line to estimate the surface area of that section. The volume of the 
section is then calculated by taking the surface area multiplied by 



the thickness of the section as it was originally cut on the freez- 
ing microtome. The density of cells or neurons is then calculated 
by dividing the number of cells or neurons estimated using the 
optical fractionator by the volume of each section determined by 
the Cavalieri method. The 10 sections are then ordered by den- 
sity. In contrast, the surface area of each flattened cortical sample 
was measured using ImageJ and multiplied by the depth of VI 
cortex measured in the intact left hemisphere cut as brain slices to 
determine the volume of each sample. The average of the isotropic 
fractionator and flow fractionator estimates was then divided by 
the volume of each sample to calculate cell or neuronal density 
(Table 1 ) . In addition, these averages were ordered by density and 
placed into 10 bins of approximately six samples each (the 10th 
bin contained seven samples). The non-parametric T statistic was 
then calculated on the ordered (a) flattened cortical sample and 
(b) intact brain section density estimates (Table 2). Pursuant to 
the third aim of determining the distribution of particles in the 
z-domain of brain slices from tissue fixed before sectioning while 
frozen on a sliding microtome, we normalized the z-axis location 
of each cell within the corresponding depth of the sampling site 
in which it was observed and divided this normalized data into 10 
bins corresponding to percentiles of tissue depth. We then used 
linear regression to investigate the presence of a statistically dis- 
tinguishable relationship between the frequency and location of 
cell counts along the z-axis. 

RESULTS 

Our investigation of chimpanzee (Pan troglodytes) primary visual 
cortex (VI, Figure 1) revealed excellent agreement between the 
estimated number of cells or neurons derived from tissue pro- 
cessed as samples from flattened cortex and quantified manually 
using the isotropic fractionator (Herculano-Houzel and Lent, 
2005) and automatically using the flow fractionator (Collins et al., 
2010; Young et al., 2012), or from tissue processed as coronal 



Table 2 | Summary statistics table. 



Analysis Tissue Quantification Parameter Test statistic P-value Sample Figures 



Within 



Between 



Flat 




Flow cytometry 


Total cells 


Nonparametric paired T-test 


0.486 


n 


= 61 










Lin's concordance statistic 


0.877 


n 


= 61 


Flat 




Neubauer chamber 


Total cells 


Nonparametric paired T-test 


0.221 


n 


= 61 


Slice 




Stereology 


Total cells 


Gundersen CE {m = 1 ) 


0.027 


n 


= 324 








Total neurons 


Gundersen CE (m = 1) 


0.018 


n 


=316 








Volume (Nissl) 


Gundersen CE {m = 1 ) 


0.017 


n 


= 324 








Volume (NeuN) 


Gundersen CE {m = 1 ) 


0.008 


n 


= 316 


Flat vs. 


Flat 


Flow cytometry vs. Neubauer chamber 


Total cells 


Nonparametric paired T-test 


0.239 


n 


= 61 










Lin's concordance statistic 


0.892 


n 


= 61 


Flat vs. 


Flat 


Flow cytometry vs. Neubauer chamber 


Percent neuron 


Nonparametric paired I-test 


0.124 


n 


= 61 










Lin's concordance statistic 


0.416 


n 


= 61 


Slice vs 


. Flat 


Stereology vs. Isotropic fractionator 


Cell density 


Nonparametric paired T-test 


0.922 


n 


= 10 










Lin's concordance statistic 


0.787 


n 


= 10 








Neuron density 


Nonparametric paired 7"-test 


0.027 


n 


= 10 










Lin's concordance statistic 


0.599 


n 


= 10 



Figure S1 
Figure S2 



Figure S3 
Figure S4 
Figure 2 
Figure 2 



Summary table showing the statistical tests performed to assess variation in measurements taken with a single quantification method (Within) and between 
quantification methods (Between) in tissue processed as flattened cortical samples (Flat) or brain slices (Slice). P-values lower than 0. 05 were considered significant. 
See Methods for full details of statistical analysis. 
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FIGURE 2 | Summary cellular density plots. Summary plots showing the 
density (in thousands of cells per mm 3 ) of cells and neurons estimated from 
flattened cortical samples (Flat) and brain slices (Slice) (A,B). In these plots, 
the flattened cortical sample estimate is the average of counts using the 
Neubauer chamber and flow cytometry for each sample. Estimates from 
flattened cortical samples were ordered and placed into 10 bins of 6 
sequential values (the 10th bin contains 7 values, n = 61 ) and are shown for 
all cells by an open triangle and for neurons by a closed triangle (A). 
Stereological density estimates from brain slices (Slice) were ordered 
(n = 10) and are shown for all cells by open circles and for neurons by closed 
circles (A). (B) Is a boxplot depicting the median (thick line), inter-quartile 
range (box), minimum (lower dotted error bar), and maximum (upper dotted 
error bar) cellular densities in flattened cortical samples and brain slices. 



sections and quantified using the computer-assisted optical frac- 
tionator method (Schmitz and Hof, 2005, Figure 2). Specifically, 
the stereological estimate (961 million) of total cell number from 
coronal slices of VI was 94.6% of the estimate based upon manual 
and automated counts (1 billion) of flattened cortex. The percent 
of all the cells that are neurons, calculated from flattened sam- 
ples (65.4%), was 90.3% of the percent estimated using stereology 
(72.4%). Manual and automated count estimates (658 million) 
of the number of neurons in VI were 95.6% of the stereological 
estimate (695 million, Table 1). 

The current investigation also found that VI volume of the 
right hemisphere estimated from images of flattened samples 
using ImageJ (5.3 cm 3 ) was within 10% of VI volume in the 
left hemisphere estimated from serial sections using the Cavalieri 
stereological method (4.8 cm 3 , Table 1). Thus, the density of 
cells in flattened samples of VI ranged from 110,783 to 329,632 
cells/mm 3 , with an average of 193,867 cells/mm 3 (Figure 2, 
Table 1). The density of neurons in flattened samples of VI 
ranged from 52,936 to 234,039 neurons/mm 3 , with an average of 
126,823 neurons/mm 3 . The density of cells in coronal slices of VI 
ranged from 123,263 to 213,806 cells/mm 3 , with an average of 
197,331 cells/mm 3 . The density of neurons in coronal slices of VI 
ranged from 108,656 to 162,511 neurons/mm 3 , with an average 
of 146,008 neurons/mm 3 . These results appear in Figure 2 and 
Table 1. Therefore, our investigation revealed that the density of 
cells and neurons from flattened samples were 98.2 and 86.9% of 
the stereological estimates, respectively. 

Corroborating the results of previous studies (Collins et al., 
2010; Young et al., 2012), we observed excellent agreement 
between repeated measurements using the isotropic fraction- 
ator (non-parametric T-test, P = 0.221, Lin's Concordance = 
0.812, n = 61) and the flow fractionator (non-parametric T-test, 
P = 0.486, Lin's Concordance = 0.877, n = 61) to estimate the 
total number of cells. Similarly, we observed low within-sample 
variance when using the optical fractionator to estimate the num- 
ber of neurons [Gundersen coefficient of error (CE) = 0.018 
(Gundersen and Jensen, 1987); 1st CE = 0.028, and 2nd CE = 
0.021 (Schmitz and Hof, 2005)] and the total number of cells 
[Gundersen CE = 0.027 (Gundersen and Jensen, 1987); 1st CE 
= 0.023, and 2nd CE = 0.018 (Schmitz and Hof, 2005)], or when 
using the Cavalieri method to estimate VI volume (Gundersen 
CE = 0.012, Gundersen and Jensen, 1987). In addition, we 
observed no difference when comparing estimates of the num- 
ber of neurons (non-parametric T-test, P = 0.124, n = 10) or the 
total number of cells (non-parametric T-test, P = 0.239, n = 61) 
based on manual counts under the microscope and automated 
counts using the flow cytometer (Table 1). 

Our investigation employed a constant estimate of the corti- 
cal depth of VI (~ 1.8 mm, n = 42) to determine the volume of 
samples from flattened cortex because estimates based upon sam- 
ple weights were more variable, likely as a result of the inclusion 
of some white matter. As noted, this approach produced a vol- 
ume estimate of right VI that was 10% larger (5.3 mm 3 ) than 
the (4.8 mm 3 ) estimate obtained using the Cavalieri method in 
brain slices from left VI. Despite this difference, our calcula- 
tions of total cellular density in chimpanzee VI comparing an 
average of the isotropic fractionator and the flow fractionator 



with the optical fractionator were statistically indistinguishable 
(non-parametric T-test, P = 0.922, n = 10). In contrast, our cal- 
culations of neuronal density using these techniques were statisti- 
cally distinguishable (non-parametric T-test, P = 0.027, n = 10, 
Table 2). The increased estimate of VI volume from flattened cor- 
tex highlights the need for the anatomical specificity available in 
stained brain slices and may account for the discrepancy between 
average neuronal density estimates using these techniques (see 
Discussion). 
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FIGURE 3 | Distribution of particles along the z-axis of Nissl-stained 
brain slices. Histogram of data based on 1399 cells (nuclei) from 207 
sampling sites throughout V1 illustrating the frequency of cell counts (Y 
axis) normalized over the depth (X axis) of sections fixed before cutting 
frozen on a sliding microtome and mounted before staining for Nissl 
substance. The convention here is that bins to the left (e.g., bin #0.1) are 
near the coverslip or the exposed edge of tissue during staining 
(superficial), and bins to the right (e.g., bin #1 ) are near the glass slide 
(deep). 



The current work revealed that the distribution of particles 
along the z-domain of sectioned materials fixed before section- 
ing while frozen on a sliding microtome and mounted on a glass 
slide before staining for Nissl-substance was vulnerable to the 
"lost caps" phenomena (Baryshnikova et al., 2006), and exhib- 
ited a linear decrease in cell density with increasing tissue depth 
(superficial = near coverslip, deep = near glass slide, Figure 3). 
Specifically, our measurements indicated that the top and bottom 
10% of tissue (bins 0.1 and 1, respectively) each exhibited ~40% 
of the average z-domain cell count ("lost caps," see Baryshnikova 
et al., 2006). Including the entire z-axis in the stereological probe 
would therefore have resulted in an underestimate of the cellu- 
lar population by ~12.9%. After excluding the top and bottom 
10% of tissue, we identified a linear decrease in the density of cell 
counts as the depth of tissue increased, such that tissue near the 
slide (deep) contained ~20% fewer counts than (superficial) bins 
near the coverslip [f ( i, 6) = 7.7, P = 0.032, Adjusted R 2 = 0.48]. 
The current work, based upon the results of a pilot study indi- 
cating the use of 2 |xm guard zones, effectively sampled the tissue 
from 2 to 12 |im of the average 16.8 u,m available, or from the 
12th to 71st percentiles along the z-axis. Based on the average 
counts per bin after excluding the top and bottom 10 percentiles 
(bin #0.1 [adjacent to the coverslip], and bin #1 [adjacent to the 
glass slide]), our calculations indicate that this sampling range 
(12-71%, or bins 2-7) resulted in an estimate of cell number 
approximately 2.6% above the sample mean. 

DISCUSSION 

Our estimates of cellular and neuronal populations in the pri- 
mary visual cortex of a chimpanzee using the isotropic frac- 
tionator, the flow fractionator or the optical fractionator agree 
within 5% (Table 1). This work also contributes data indicating 
that the measurement of VI volume from flattened preparations 
agrees with stereological estimates from brain slices within 10% 
(Table 1). Our investigation therefore provides evidence of the 



agreement between estimates of the density of cells and neurons 
based upon sampling flattened cortex or brain slices (Figure 2). 

We provide data on the volume of chimpanzee VI 
(4.8-5.3 cm 3 ) consistent with the 4.64 cm 3 (De Sousa et al., 
2010) and 5.52 cm 3 (Bush and Allman, 2004) estimates from pre- 
vious studies. In particular, our stereological estimate of 4.8 cm 3 
supports the De Sousa et al. (2010) study that also used the 
Cavalieri method to estimate the volume of VI in the left hemi- 
sphere of chimpanzees at 4.64 cm 3 . The difference between our 
estimates in the left (4.8 cm 3 ) and right hemispheres (5.3 cm 3 ) 
may reflect inter-hemispheric asymmetry in VI volume, given 
that previous work has reported differences in humans of 5% 
(Murphy, 1985) and 13% (Andrews et al, 1997). This difference 
in VI volume may be related to behavioral specializations as 
has been observed in the primary motor cortex related to hand 
preference (Sherwood et al., 2007a,b; Hopkins et al, 2010). On 
the other hand, it is possible that mechanical deformation of the 
cortical sheet induced during the flattening procedure increased 
our surface area measurement which, when combined with depth 
measurements from the intact left hemisphere, produced a slight 
over-estimate of VI volume in the right hemisphere. However, 
given that our estimates were within 10% of each other, and that 
our calculations fall well within the range suggested by previous 
studies and exhibit inter-hemispheric asymmetry similar to what 
has been reported in humans, these results suggest that our 
measurements provide reasonable estimates of the volume of VI 
in the chimpanzee brain. 

Results from our investigation of the distribution of parti- 
cles along the z-axis in sectioned materials stained for Nissl after 
mounting suggest that the use of guard zones to exclude the top 
and bottom 10% of tissue avoids underestimating the total num- 
ber of cells by approximately 13% (effect of "lost caps," Figure 3, 
see Baryshnikova et al., 2006). However, the implementation of 
2 (Am guard zones resulted in an effective probe sampling range 
of the 12— 71st percentiles along the z-axis which resulted in a 
slight overestimation of cell number by approximately 2-3%, 
based upon the observed decrease in cell density as the observer 
approaches the glass slide. However, previous work has reported 
an increase in cell density as the observer moves from the tissue 
near the coverslip toward the glass slide upon which the tissue 
was mounted, and the absence an effect of "lost caps," such that 
cells near sectioned edges were not notably fewer in number than 
in the middle of the brain slice (Carlo and Stevens, 2011). The 
inverse effects observed in these two studies are likely explained 
by the observation that the tissue Carlo and Stevens (201 1) mea- 
sured was processed in a manner identical to that employed in the 
current investigation, except for two things: (1) Carlo and Stevens 
(2011) report flipping the tissue such that the edge of the tissue 
immediately adjacent to but above the knife blade during sec- 
tioning was positioned closest to the coverslip in the final stained 
product, whereas the current study did not flip the tissue in this 
manner, and (2) Carlo and Stevens (2011) measured brain tis- 
sue that underwent perfusion fixation prior to sectioning, while 
the current investigation measured tissue that underwent a saline 
flush during perfusion and was subsequently fixed by immer- 
sion before sectioning. The opposing gradients of cell loss are 
therefore likely the result of the same phenomena in both studies 
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such that cells immediately above the microtome blade during 
sectioning exhibit an increased density, an effect of tissue com- 
pression that occurs prior to mounting the tissue (Baryshnikova 
et al., 2006), and are simply being measured across the tissue in 
opposing directions. Similarly, the presence of "lost caps" in the 
current investigation but the absence of this effect in Carlo and 
Stevens (2011) data is therefore likely the result of whether the 
brain was fixed by perfusion (resulting in no lost caps), or not 
(resulting in approximately 60% cell loss at the cut edges). This 
conclusion is supported by similar observations from different 
labs that cells are lost at cut edges in immersion fixed (Andersen 
and Gundersen, 1999; Bahney and von Bartheld, 2014) but not 
perfusion fixed tissues (Gardella et al, 2003). Thus, these stud- 
ies together suggest that compression occurs regardless of fixation 
and is characterized by an effect somewhere between the currently 
reported 20% (clearance angle = 19°) and previously reported 
50% (clearance angle = 7.5°, Carlo and Stevens, 2011) decrease 
in cell density as the observer moves away from the edge of tis- 
sue adjacent to and above the knife edge during sectioning. In 
addition, perfusion fixation appears to mitigate the effect of "lost 
caps," which results in the loss of a significant number of cells or 
particles near the surface of sections from immersion-fixed tissues 
(~60%). Given that the ultimate goal of stereological investiga- 
tions is to estimate biological parameters of interest in such a way 
that the original relationships among those parameters are not 
systematically distorted, our data suggest that the use of guard 
zones in these tissue preparations are critical, and that awareness 
of the extent to which the stereological probe employed samples a 
known distribution of objects within the available tissue can pro- 
vide reasonable estimates, for example, of the number of cells and 
neurons. In summary, our data provide evidence in support of 
the traditional stereological approach to focus on sampling the 
"core" of sectioned material, because this avoids the largest source 
of bias due to the "lost caps" effect at tissue edges (von Bartheld, 
2012). In addition, our data provide evidence to suggest that the 
use of guard zones resulting in an effective sampling range from 
anywhere within the 10-30th to the 70-90th percentiles produce 
results within 3-4% of the observed or sample mean. 

This study corroborated previous work (Young et al, 2012) 
comparing the total number of cells in a sample measured using 
the isotropic fractionator and the flow fractionator (Figures Sl- 
S3). We found that the total number of cells estimated from 
flattened cortical samples was 5% greater than the stereological 
estimate (Table 1). However, our stereological investigation likely 
overestimated the actual number of cells in our sample by approx- 
imately 3%, resulting in a discrepancy between hemispheres of 
approximately 8%. Given that previous investigations of human 
VI, as well as the current volumetric analysis, indicate inter- 
hemispheric differences on the order of 5-13% (Murphy, 1985; 
Andrews et al, 1997), these data suggest that the observed dif- 
ferences in cell number reflect biological asymmetry between the 
hemispheres. The close agreement between these results and pre- 
vious studies therefore strongly suggests that each of these distinct 
techniques estimate a similar parameter, the total number of cells 
in chimpanzee VI. 

This work provides evidence in agreement with previous stud- 
ies investigating the percentage of neuronal cells in samples 



measured using the isotropic fractionator and the flow fraction- 
ator (Collins et al, 2010; or Figure S4). It should be noted that 
the estimate of the absolute number of neurons in tissue samples 
from flattened cortex is a derivative of the estimate of the total 
number of cells using DAPI. This contrasts with the stereological 
estimate of neuron number, which is based upon an independent 
run of the optical fractionator in a series of sections stained for 
NeuN, adjacent to the Nissl-stained series. The percentage of neu- 
ronal cells based upon our stereological investigation is calculated 
as the number of cells in NeuN-stained sections divided by the 
number of cells in Nissl-stained sections. Therefore, our compar- 
ison of the number of neurons using these techniques is two-fold: 

(a) we compared the percentage of all the neuronal cells, and 

(b) we compared the absolute number of neurons estimated. We 
found that the mean percentage of neuronal cells in chimpanzee 
VI estimated from flattened cortical samples was 65.4%, while 
the mean percentage of neuronal cells estimated using stereology 
was 72.4%. Given that our stereological analysis likely produced 
a ~3% overestimate of cell number, the percentage of neuronal 
cells estimated from brain slices may be approximately ~74%. 
This difference in the percentage of neuronal cells is likely due 
to the inclusion of a minimal amount of white matter in the 
flattened cortical sample preparation. Next, we calculated the 
absolute number of neurons and discovered that the estimate 
from flattened cortical samples (658 million) was 5% lower than 
the stereological estimate of 695 million (Table 1). Given that pre- 
vious evidence indicates a >5% rightward bias in VI size and cell 
number, the observation that right VI contains 5% fewer neu- 
rons suggests that neurons may be lost during implementation 
of the isotropic fractionator method (Clarke and Oppenheim, 
1995). Although future research may further clarify this issue, 
the reasonably close (~10%) agreement between these estimates 
suggest that our quantification procedures approximate the same 
parameter, the total number of neurons in chimpanzee VI. 

Our data provide evidence that there are approximately 
195,000 (193,867-197,331) cells per mm 3 of chimpanzee VI, out 
of which around 136,000 (126,823-146,008) are neurons. These 
data differ somewhat from the results of a previous study of cellu- 
lar density in chimpanzee VI which reported ~330,000 cells per 
mm 3 , out of which 208,930 were neurons (Lewitus et al., 2012). 
Interestingly, both of these results place the density of neurons in 
chimpanzee VI in range of the upper 60 percentiles. However, 
the discrepancy between the absolute density of neuronal and 
non-neuronal cells in these studies is likely due to a combination 
of factors including tissue preparation, the type of fixation and 
the histological staining protocol employed. For instance, Lewitus 
et al. (2012) excluded layer 1 from their estimates as a result 
of variable tissue quality, and measured cortical thickness only 
twice for each section investigated, the combination of which may 
have contributed to the higher cellular densities they observed. 
However, our results are in agreement with this and other pub- 
lished studies indicating a high density of neurons in primary 
visual cortex, particularly in primates (Collins et al., 2010, 2013; 
Campi et al., 2011; Lewitus et al., 2012). 

The density of cells and neurons based on flattened samples 
were 98.2 and 86.9% of the stereological estimates, respectively. 
However, the observation that our stereological data probably 
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overestimated (~3%) the total number of cells may account for 
the observed ~2% difference in cell density between hemispheres. 
These results demonstrate that the density of cells can be deter- 
mined with high precision, whether measured using traditional 
stereological tools or with more recently developed techniques 
such as the isotropic fractionator and flow fractionator. In our 
hands, these distinct methodologies produced estimates of the 
density of neurons that agreed within 15%. The source of this 
discrepancy may involve the inclusion of a minimal amount of 
white matter in flattened cortical samples during dissection, the 
slight overestimation of regional volume because of deformations 
in the cortical sheet induced during the flattening procedure, the 
loss of neuronal cells, or some combination of these. Despite 
these relatively small differences, our ability to replicate the com- 
position and density of cellular populations in the chimpanzee 
primary visual cortex provides evidence supporting our ability to 
compare results from distinct methodologies in determining the 
distribution of different types of cells in the brain. 

CONCLUSIONS AND SIGNIFICANCE 

Our research indicates the reproducibility of estimates of the cel- 
lular composition of brain tissue using distinct tissue processing 
and quantification methodologies when the accurate identifica- 
tion of a single cortical region of interest is possible. Our results 
thus demonstrate the validity of using the isotropic fractionator, 
flow fractionator, or optical fractionator to obtain cell and neuron 
counts. Ultimately, these results facilitate a more comprehensive 
understanding of neurobiological circuitry by presenting inves- 
tigators with evidence that these techniques estimate the same 
parameters of interest, allowing the investigator an additional line 
of evidence to aid in the interpretation of experimental results. 

Our results are also significant because they indicate that VI in 
chimpanzees, our closest living genetic relatives, conforms to the 
general finding in primates that VI contains more densely packed 
neurons than other cortical areas (Rockel et al., 1980; Carlo and 
Stevens, 2013; Collins et al., 2013; Young et al., 2013). This means 
that, beyond the observation that layer 4 contains small neu- 
rons, each of the cortical layers in VI contains neurons that are 
on average smaller than in other cortical areas. Smaller neurons 
tend to have fewer synaptic inputs (Jacobs et al., 2001; Elston 
et al., 2005) and are thus thought to be less likely to dramatically 
transform incoming sensory signals. Primate brains are known 
to contain a large proportion of neocortical tissue devoted to 
visual perception. Furthermore, compared to other mammals, VI 
in primates has the specialized task of distributing sensory input 
to cortical regions throughout the visual system that are integral 
for visual percepts. Accordingly, our findings are consistent with 
the observation that primate VI is characterized by small, densely 
packed neurons, which may have facilitated the evolution of pri- 
mate VI into a region particularly effective at providing multiple 
downstream cortical regions with the relatively raw sensory input 
necessary for higher-order functional operations. 
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Figure SI | Analysis of flow cytometry estimates. Repeated measures 
graph depicting the agreement between the first (red) and second (blue) 
density measurements on ordered samples (A). The vertical dotted lines 
indicate every 5th sample (A). Bland-Altman plot showing the agreement 
between the first and second measurement taken from each sample by 
comparing the difference between repeated measurements on the Y axis 
against the mean of each sample estimate plotted along the X axis (B). 
The blue dotted line indicates the sample mean, the green dotted lines 
indicate the 95% interval, and the red dotted line indicates a single 
standard deviation from the mean (B). All density estimates are in 
thousands of nuclei per mm 3 (A,B). 

Figure S2 | Analysis of manual counts using the Neubauer chamber. 

Repeated measures graph depicting the agreement between the first 
(red) and second (blue) density measurements on ordered samples (A). 
The vertical dotted lines indicate every 5th sample (A). Bland-Altman plot 
showing the agreement between the first and second measurement 
taken from each sample by comparing the difference between repeated 
measurements on the Y axis against the mean of each sample estimate 
plotted along the X axis (B). The blue dotted line indicates the sample 
mean, the green dotted lines indicate the 95% interval, and the red dotted 
line indicates a single standard deviation from the mean (B). All density 
estimates are in thousands of nuclei per mm 3 (A,B). 

Figure S3 | Analysis of flattened cortical sample estimates. Repeated 
measures graph depicting the agreement between the average density 
count using the Neubauer chamber (red) and the flow cytometer (blue) on 
ordered samples (A). The vertical dotted lines indicate every 5th sample 
(A). Bland-Altman plot showing the agreement between these 
quantification procedures by comparing the difference between estimates 
on the Y axis against the mean of each sample estimate plotted along the 
X axis (B). The blue dotted line indicates the sample mean, the green 
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dotted lines indicate the 95% interval, and the red dotted line indicates a 
single standard deviation from the mean (B). All density estimates are in 
thousands of nuclei per mm 3 (A,B). 

Figure S4 | Analysis of neuronal percentage. Repeated measures graph 
depicting the agreement between the average neuronal percentage using 
the Neubauer chamber (red) and the flow cytometer (blue) on ordered 
samples (A). The vertical dotted lines indicate individual samples (A). 
Bland-Altman plot showing the agreement between these quantification 
procedures by comparing the difference between estimates on the Y axis 
against the mean of each sample estimate plotted along the X axis (B). 
The blue dotted line indicates the sample mean, the green dotted lines 
indicate the 95% interval, and the red dotted line indicates a single 
standard deviation from the mean (B). All density estimates are in 
thousands of nuclei per mm 3 (A,B). 
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